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Abstract 

We prove an area law with a logarithmic correction for the mutual information for fixed points of local 
dissipative quantum system satisfying a rapid mixing condition, under either of the following assumptions: 
the fixed point is pure, or the system is frustration free. 


1 Introduction 

One of the problems common to quantum information and condensed matter physics is understanding 
correlations and entanglement in many-body quantum states. This is motivated by the observation that 
many interesting states behave very differently from random states: while in the latter the entanglement 
entropy of a sub-region scales as the volume of the region, in the states that occur in quantum many-body 
systems this entanglement entropy is often seen to scale only with the boundary of the region. This surprising 
behaviour has been called the area law (the terminology comes from 3D systems), and it is a well-studied 
conjecture that ground-states of (gapped) local Hamiltonians should satisfy an area law. Even in the case 
of critical systems and gapless Hamiltonians, evidence suggests that the groundstate still has a sub-volume 
growth of the entanglement, with a rate proportional to the boundary of the region times a logarithmic 
correction [30, 38] * 1 . While it has been formally proven only for ID systems [2, 3, 7, 8, 17], area laws have 
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1 Some authors call such situation a violation of the area law, while others speak of an area law with a logarithmic correction. 
We will use the second terminology, but we will use the term sub-volume law for any faster growth rate of entanglement that is 
slower than the volume of the region. 
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also been proven in specific cases in higher dimensions (harmonic lattice systems [35], models satisfying Local 
Topological Quantum Order [32], perturbations of gapped Hamiltonians satisfying an area law [28, 31], and 
with a logarithmic correction for fermionic systems [12, 43]), and are the subject of active research. 

Recently, a different class of states arising in quantum many-body systems has been attracting attention 
in the quantum information literature: fixed points of (local) dissipative processes. More precisely, fixed 
points of semigroups of trace preserving, completely positive linear maps. The motivation is two-fold: on the 
one hand, such processes model most of the different types of noise that can be found in nature, and therefore 
provide a more realistic model for physical systems, since in practice no system will be completely isolated. On 
the other hand, proposals have been made to artificially engineer such dissipative interactions in order to have 
a determined quantum state as a fixed point, effectively making them “dissipative machines” for producing 
useful/interesting quantum states [23, 41]. This dissipative state engineering has been experimentally shown 
to be a robust mechanism to maintain coherence [5, 24]. 

A natural question then arises: is there an area law (either strict or with a logarithmic correction) in this 
context? Note that, since fixed points of dissipative evolutions are generically not pure, we must find another 
measure of entanglement or correlations, since local entropy is no longer a useful measure for mixed states (as 
the trivial example of the maximally mixed state shows). [44] proposed instead using the mutual information, 
a measure of correlations between two parts of a quantum state. It has the advantage of coinciding with 
entanglement entropy for pure states, and it upper bounds operational measures of entanglement in the 
mixed-state case, such as the distillable entanglement [6]. In [44] it was also shown that thermal states of 
local Hamiltonians satisfy an area law for the mutual information. However, since thermal states do not cover 
all the possible fixed points of dissipative systems. Pinning down under which conditions an area law for the 
mutual information holds for general dissipative quantum many-body systems is thus an interesting open 
problem. 

For Hamiltonian systems, the main assumption that is usually made is the presence of a spectral gap: a 
non-vanishing separation between the two lowest energy levels of the Hamiltonian. In the dissipative setting, 
instead of spectral assumptions, it is more natural to make assumptions on the speed of convergence of 
the dissipation towards its fixed point (a quantity that is not controlled by the spectrum alone [39]), or 
equivalently on the so-called mixing time. In this work we restrict to systems for which the mixing time 
scales logarithmically with the system size. In a previous paper [10] some of the authors showed that such 
systems - which we called rapid mixing - are stable under local perturbations. Similarly to the gap of closed 
quantum systems, proving rapid mixing for dissipative systems is a daunting task. There are however some 
interesting key examples: state preparation of graph states [21], classical Glauber dynamics for the Ising 
model in 2D (in some range of parameters) [27, 29], among others. 

Rapid mixing is implied in many cases by a well studied property of dissipative evolutions, namely the 
existence of a system-size independent Log-Sobolev constant for primitive reversible Liouvillians [14-16, 22, 
45]. Under such assumptions, in [20] it was proven a bound on the mutual information of the form: 

I(A : A c ) < cloglog ||p _1 || \dA\. (1) 

In the latter bound p is the fixed point of the evolution, and therefore p _1 will usually depend at least 
exponentially on the total system size, and sometimes even worse. As recognized by the authors of [20], this 
poses a serious problem in considering Eq. (1) a satisfactory area law. It indicates however that rapid mixing 
seems to be the required condition to have an area law in the dissipative setting. This is exactly what we 
prove in the present paper, with the following results: 

1. if the system satisfies rapid mixing and the fixed point is pure, then it satisfies an area law with a 
logarithmic correciton for the entanglement entropy; 

2. if the system satisfies rapid mixing and is frustration free, meaning that the local terms of the Liouvillian 
share a common steady state, then such fixed point satisfies an area law with a logarithmic correciton 
for the mutual information. 

Compared with (1), the bounds we obtain do not have any dependence on the total system size. Moreover 
we do not require primitivity or reversibility of the generators of the evolution, and we only require rapid 
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mixing instead of a system-size independent Log-Sobolev constant (a strictly weaker assumption, since the 
Log-Sobolev constant is undefined for non-primitive Liouvillians). 

It is known that there is a connection between area laws and decay of correlations [7, 44]. Therefore it 
does not come as a surprise that with the tools we have developed for proving the area law we can also prove 
a decay of correlations measured with the mutual information. It is worth noting that, with the results 
available in the literature and to the best of our knowledge, it is not possible to derive the area law from the 
type of decay of correlations we will show here. 

As we have mentioned earlier, for groundstates of closed systems a logarithmic correction is usually 
considered a signature of a gapless Hamiltonian. For open systems and mixed states the situation is less 
clear: already in [44] it was shown that thermal states of local Hamiltonian satisfy an area law without a 
logarithmic correction irrespective of the gap of the Hamiltonian. For these states the bound we obtained is 
therefore not optimal. We do not know whether there exists systems which saturate our bound, or if instead 
the correction is only an artefact of the proof. If there exist systems that saturate our bound, it would then 
imply that their fixed point have a very interesting property: while still satisfying an exponential decay of 
correlations, they do not satisfy an area law without a logarithmic correction. Having an example of such 
state which can also be efficiently prepared with a dissipative process would be interesting on its own, as it 
could lead to new insight on the relationship between area laws and decay of correlations. 

We conjecture that rapid mixing alone, without any additional assumptions, should imply an area law for 
mutual information, but we do not have a formal proof. The fact that we have two different proofs of an 
area law, requiring different extra assumptions (pure fixed point on the one hand, frustration freeness on the 
other) is strong evidence for this conjecture. 

The paper is organized as follows. In section 2 we set up the problem and introduce the necessary notation 
and definitions. In section 3 we prove two lemmas regarding localization properties of the fixed point of the 
dissipative maps: Lemma 11 is based only on the rapid mixing assumption, and will be used in section 4 to 
prove decay of correlations and in section 5 to prove the area law for pure fixed points; Lemma 12 instead 
requires the extra assumption of frustration freeness, and will be used in section 5 to prove the area law for 
the mutual information in the case of mixed fixed points. 


2 Setup and notation 

Let Bab = Ba ® 'Hb be a finite dimensional complex Hilbert space, representing a bipartite quantum system. 
We denote by B{Bab) the space of bounded linear operators on Bab- A state is given by a positive semi- 
definite operator pab 6 B(Bab), normalized to have trace equal to one. The reduced density matrix of the 
subsystem A (resp. B) will be denoted by pa (resp. ps), and it is given by pa = Tr^ pab (resp. ps = Tr^ pab), 
where the partial trace Tr^ is defined to be the unique linear operator Tr^ : B(Bab) -*■ B(Bb) such that 
Tr^rr® y) = yTi(x) for all x in B(Ba) and all y in B(Bb) (Tr s is similarly defined). 

We will use the standard Dirac notation for Hilbert spaces, denoting vectors as \<f>), adjoint vectors as {4\, 
((f>\il>) for the scalar product, and |0)('01 for rank one linear maps. The canonical basis will be indexed by 
natural numbers starting from zero: |0), |1),... |n). 

We will denote by poly(x) any polynomial in the variable x with real coefficients and arbitrary degree. 

2.1 Measures of correlations 

Given a state pab 6 B(Bab) of a bipartite system, there are a number of possible measures of how “distant” 
the state pab is from being a product state, i.e. of the form x ® y for some x in B(Ba) and y in B(Bb)- 
Since a product state represents a system in which measurements over the subsystem A are independent of 
measurements over the subsystem B , we will talk of correlation measures between subsystem A and B. We 
will need to define and use three of such measures. We will follow the same terminology of [20]. 

Definition 1 (Correlation measures). 
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Covariance correlation ; 


C(A:B)= max \(M ® N) - (M) (N)\ 
MzB(Ha),N£B('Hb) 

||M||<1,|| AT||<1 


max \Tr [M ® N(pab ~ PA® Pb)]\; 
M£B(Ha),N£B(Hb ) 

]|M||<l,||iV||<l 


where (O) = Tt(Opab) is the expectation value of the observable O acting on pab- 
• Trace distance correlation: 


T(A:B)= max |Tr [F{p AB - p A ® Pb)]\ 
FsBCHab) 

\\F\\<1 

= |PAB ~ PA ® Pb\i ■ 


• Mutual information correlation: 


I(A: B) = S(p A ) + S( PB ) - S(pab); 

where S(p) = - Tr(plog 2 p) is the von Neumann entropy of the state p. 

When it is not be clear from context which state pab we are considering, we indicate it in a subscript to 
avoid ambiguity, and write C(A : B ) p , T{A : B) p and I(A : B) p . 

As should be clear from the definition, C(A : B ) is always upper bounded by T(A : B). Moreover, by 
Pinsker’s inequality [34]: 

C(A : B) < T(A : B) < 2 y/l(A:B). 

Therefore, mutual information is the strongest correlation measure. It is also a well known consequence of 
the Alicki-Fannes-Audenaert inequalities [1, 4, 11] that there is a non-linear inverse relationship between 
trace distance and mutual information: there is a differentiable function f(x) vanishing at zero such that 
I (A: B) < f(T(A : B)). Such non-linear equivalence between the two measures will allow us to take bounds 
on T(A : B ) (which in the context of our assumptions will be easier to deduce) and obtain information on 
the behavior of I(A: B). 

We state this result in a form that will be more convenient for us and, for the sake of completeness, 
present a short proof of it. 

For x e [0,1], hb(x) = -aTog 2 x - (1 - x) log 2 (l - x) denotes the binary entropy function. 

Theorem 2. The following inequalities hold: 

Fannes-Audenaert [4, 11] 

Let p,a e B(€ d ), and let S = \p - cr^ < 1. Then 

|S(p) - 5(a)] < 2<51og 2 (d- 1) + 2 h b (S), (2) 

Alicki-Fannes [1] 

Let p AB ,a AB eB(€ 4 ®C dfl ), and let 6 = \\p AB - cr AB \\ 1 < 1. Then 

| S(p AB \p B ) - S(a AB \a B )\ < 4(5log 2 d A + 2 h b (S), (3) 

where the conditional Von Neumann entropy S(p AB \p B ) is defined as S(p AB \p B ) = S(p AB ) - S(p B ). 
Combining the two previous inequalities, we will obtain the desired non-linear bound on I(A : B). 
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Corollary 3. Let 5 = \\p AB ~ & AB < 1. Then 

\I(A : B) p - I {A : B) a | < 65 log 2 d A + 4 h b (S). 

In particular, if we take <j ab = p A ® p B , we have that I(A : B) a = 0,5 = T(A : B) p , and thus 

I {A : B) p < 6T(A : H) p log 2 d A + 4h b (T(A : H) p ). (5) 

Proof. Note that |/cA B - <t‘ 4b || 1 < 5 => \\p A - cr' 4 | < 5. Applying triangle inequality and equations (2) and (3) 
gives 

I I(A : B) p - I(A : B ) a I < |S(p A ) - S(rr A )| + \S(p AB \p B ) - S(a AB \a B )\. 

□ 


2.2 Many-body quantum systems 

Let us now recall the standard definitions and the common notation for many-body quantum systems. 

We will consider a quantum system defined on the square lattice T = iP equipped with the graph metric, 
where at each site x e T we associate a finite-dimensional complex Hilbert space 7i x . We choose to work with 
a square lattice for simplicity of exposition, but the results presented can be generalized straightforwardly 
to graphs with polynomial growth (i.e. with balls size growing polynomially with the diameter). The ball 
centered at x of radius r will be denoted by b x (r). We will use the following convention: given a subset A c A, 
we will denote by A(s) the smallest disjoint union of balls containing {x e A|dist(x, A) < s). For each finite 
subset A c T of the lattice, we will associate a Hilbert space 7i\ = ® X sa7L x and an algebra of observables 
Aa = 13(71 a)- We will equip Aa with the Hilbert-Schmidt scalar product ( A,B) = Tt(A*B). 

A linear map T : Aa -*■ Aa will be called a superoperator [42] to stress the fact that it is an operator acting 
on operators. Its support is defined to be the minimal set A' £ A such that T = T' <8> 1, where T' e B(Aa')- 
Positivity is defined as usual for linear maps: T is said to be positive if it maps positive operators to positive 
operators. T is called completely positive if T <8> 1 : A a ® M n -> A a ® M n is positive for all n > 1. Finally, we 
say that T is trace preserving if Tr T(p) = Tr p for all p e Aa- 

A dissipative evolution for a quantum system is given by a one-parameter continuous semigroup of 
completely positive and trace preserving superoperators {T t ■ Aa Aa}*- If p e B(H a) is the state of the 
system at time zero, then the evolution of p at time t > 0 is given by p(t) - T t (p). The assumptions on 
T t guarantee that p(t) is again a state, i.e. a positive and trace one operator. This is usually called the 
Schrodinger picture. 

We will make use of the following norm for superoperators: 

„ ,, „ „ IIT <8> l n (X)||-, 

ll T llo = sup||T® l„||i^i =sup sup - — -. (6) 

n n .Y6j4 a 8M„ IW II 1 

x*o 

Dissipative maps are contractive with respect to such norm, in the sense that ||7~| 0 < 1. 

Given a semigroup of dissipative maps {T t }, it has a generator C : Aa -*■ Aa which satisfies f^T t (p) = 
CT t (p). Such superoperator is called a Lindbladian or Liouvillian (we will use the former). The assumptions 
made on T t force a particular structure on £, which is called the Lindblad form [13, 26]. A superoperator C is 
said to be in the Lindblad form if it can be written as 

C(p) = i[p,H] + Y,LjpL* - ]-{L*Lj,p }, (7) 

j z 

where H is an Hermitian operator, (Lj)j are arbitrary operators (called the Lindblad operators), [-,-] denotes 
the commutator and the anti-commutator. We refer to [9, 42] for details on the theory of Lindblad 
operators. 

As shown in [42], given a semigroup of dissipative maps {T t }t we can define a new map representing 
the “infinite time limit” of the evolution, or in other words the projector onto the space of fixed points of 
the evolution. Tx, is again a completely positive, trace preserving superoperator, since it can be obtained as 
liniA^oo A Y,n=0 T n- 
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2.3 Uniform families 


Given a generator C, we can decompose it as a sum of local terms , i.e. terms which are still of the Lindblad 
form but with controlled support: 

C = Hz, supp£z = Z, Cz is Lindbladian. 

Zc A 

When ||>Cz|| 0 is decaying with diamZ, we will generically say that the evolution is local. More stringent 
assumptions on the decay rate of the norms of the local generators will be required, and are formalized in 
Assumptions (A-l) and (A-2) in this section. 

Since we are interested in dissipative evolutions defined on increasing sequences of lattices, and how their 
properties depend on the lattice size (often referred to as the system size), we need to define a meaningful 
way of growing the evolutions with the lattice size, adding and modifying the necessary generator terms 
appropriately. As presented in [10], the following definition of uniform families of dissipative evolutions is 
one solution to this, which is general enough to cover a wide range of models and situations. 

Definition 4. Given Acf, a boundary condition for A is given by a Lindbladian B C>A = Y,d> l B 9A , where 
supp B 9A c d d A ■■= {x e A | dist(:c, A c ) < ci). 

Definition 5. A uniform family of Lindbladians is given by the following: 

(i) infinite Lindbladian: a local Lindbladian M. defined all of TLP: Ai = T)zcz° -Mz! 

(ii) boundary conditions: a family of boundary conditions {B dA }where A = b u (L), for each u e h D and 
L> 0. 

Definition 6. A local Lindbladian C = ffz<zZ D £>z is said to be translationally invariant if Cz+ U - £ z , V« e 

7 . n . 

We say that a uniform family C = {At,B} is translationally invariant if M. is translationally invariant, 
and moreover is independent of u. 

Given a uniform family C = {At, B}, we fix the following notation for evolutions defined on A = b u (L) c F: 

C A - ^ Mz “open boundary” evolution; (8) 

Zc A 

C A = C A + C dA “closed boundary” evolution, (9) 

with the respective evolutions T A = exp (tC A ) and T A = exp (tC A ). 

Until now we have made no specific assumption on the decay rate of the norms of the local generators. 
As mentioned above, in order to meaningfully talk about locality of the evolution, we need to impose that 
\Cz\ 0 is decaying with diamZ. The rate at which such function decays classifies the system into one of the 
specific cases more usually considered in the literature: compactly supported (usually called finite-range 
interactions), exponentially decaying, super-polynomially decaying, power-law decaying, etc. We will take 
a more general approach, and will simply assume from now on that our family of Liouvillians satisfies the 
following assumptions. 

Definition 7 (Lieb-Robinson Assumptions). There exists an increasing function v(r) satisfying v(x + y) < 
v(x)v(y), such that the following conditions hold: 

sup£ llA^lo |Z|iz(diamZ) < v < oo, (A-l) 

z?x 

supsupz'(r) ^ 11^93(2:,^) II <poly(j\r). (A-2) 

x«r r d=r N 1,0 
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If I-Mzlo is exponentially decaying or is compactly supported, then one can take v(r) = exp(/ir) for some 
positive p. On the other hand, if ||A1z||^ decays only polynomially, then we must take v(r) = (1 + r) M . In the 
latter case, the Lieb-Robinson bounds, only hold if /i is bigger than a constant depending on the geometrical 
dimension of the lattice F, i.e. D if T = h D . The details of when we can apply Lieb-Robinson bounds in this 
case can be found in our previous work [10]. From now on we will simply assume that i'(-) decays sufficiently 
fast for the Lieb-Robinson bounds to apply. 

2.4 Frustration freeness 

The following definition is inspired by the analogous concept defined for closed systems and Hamiltonian 
dynamics. It captures the idea that a fixed point of a local evolution might or might not be locally steady. 
The local dissipative terms could in general have a non-trivial action on the fixed point, so that it is only the 
sum of such local effects that adds up to zero and leaves the state invariant. Assuming that the system is 
frustration free means excluding such cases. 

Definition 8. We say that a uniform family £ - {A 4,13} satisfies frustration freeness (or is frustration free,) 
if for all A and all fixed points of T t A 

M z (p°°) = 0 VZcA. (10) 

Observation 9. If £ = {M. , B} is frustration free, then each fixed point of T A is also a fixed point of T A . 

While not true in general, the frustration freeness condition is satisfied by a large class of interesting 
dissipative systems, as the following examples show: 

1. Dissipative state engineering procedures defined in [23, 41]; 

2. Locally reversible classical Markov chains; 

3. Locally detailed balanced quantum Markov processes, and in particular Gibbs samplers for commuting 
Hamiltonians [19]. 

3 Localization results 

A well known property of many-body systems with local interactions, either dissipative or Hamiltonian, is the 
existence of a finite speed of propagation. This describes how the support of a localized observable spreads in 
time during the evolution: up to an exponentially small correction, the support spreads linearly with time. 
The finite velocity at which such linear growth occurs is often called a Lieb-Robinson velocity, or sometimes 
a group velocity. (It is a property of the model and not a consequence of some relativistic effect - we are 
considering only non-relativistic models here.) While the original work focused on Hamiltonian systems and 
groups of automorphisms [25, 37], the existence of such finite speed of propagation in the lattice has been 
generalized to dissipative evolutions [33, 36], and in [10] we showed that the definition of boundary condition 
we have given allows us to recover the same type of localization properties of the evolution. 

Nonetheless, all the Lieb-Robinson localization bounds have a time-dependency, becoming worse as time 
increases, until the bound they provide becomes trivial and does not give any information at all about the 
properties of the fixed point. In the next section we want to produce results which might be interpreted 
as “infinite time” versions of Lieb-Robinson bounds. To do so, we will need to make an extra assumption 
on the evolution: we will assume that the convergence to the fixed point is fast, in the sense that scales 
logarithmically with the system size. This is formalized in the definition of rapid mixing below. Such a 
definition can be in some cases relaxed to allow convergence which is only scaling sub-linearly with respect to 
system size. (We will not pursue such generalizations here, and instead refer to [10] for guidance on which 
changes are necessary to the results below.) 

The localization lemma that we prove, Lemma 11, will be only be sufficient to prove an area law for the 
pure fixed point case, not for the mutual information. Therefore, we also prove a stronger result, Lemma 12, 
for which we will need to add the extra hypothesis of frustration freeness. 
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3.1 Rapid mixing 

In this section, we want to briefly recall a result proven in [10]. We start by recalling the definition of rapid 
mixing. 

Definition 10 (Rapid mixing). Let {Tfi } A be a family of dissipative maps, we say it satisfies rapid mixing 
if there exist c, 7 , S > 0 such that 

sup ||T, A ( ( 0 )-T^(p)|| <c|A | 5 e -‘ 7 . (11) 

p> 0 

Tr p = 1 

In analogy to the spectral gap for Hamiltonians, proving that a family of Lindbladians is rapid mixing is 
not an easy task. Nonetheless, there exists a large class of interesting systems for which we already have 
mixing time estimates that imply rapid mixing: 

1. (Trivially) non-interacting particle systems. 

2. Dissipative state engineering for graph states [21]. 

3. Quantum and classical Markov processes satisfying a Log-Sobolev inequality [22]. This includes in 
particular Glauber dynamics for the Ising model in 2D, either above the critical temperature or with 
non-zero magnetic field [29]. 

Lemma 11. Let C = {A i,B} be a uniform family of dissipative evolutions that satisfies rapid mixing, and 
suppose each T t A has a unique fixed point and no other periodic points. Fix a A and let p^ be the unique 

fixed point of T t A . Given A c A, for each s > 0 denote by p ^ the unique fixed point of T^ s \ 

Then we have: 

|Tr^(p 00 -p^)|| 1 <|A| i Ao( S ), ( 12 ) 

for some fast-decaying function Ao(s) and some positive constant S. 

The decay rate of A 0 (s) is in the same class as zz _ 1 (s), where u(s) is defined by Assumptions (A-l) 
and (A-2): it is exponential if zz - 1 (s) is exponential, polynomial if v~ 1 (s) is polynomial. In the latter case 
the degree of the polynomial controlling the decay is smaller than that of ^ _ 1 (s), but the loss is independent 
of the system size - again, this corresponds to requiring a sufficiently fast-decaying z/ - 1 (s). 

3.2 Localizing with frustration freeness 

In the following section, we want to show a property of the fixed points of a uniform family of Lindbladians 
verifying frustration freeness. We want to study the behavior of a system when it is prepared and started 
in a state, which is the fixed point of the same family but of a slightly smaller region. A reasonable guess 
would be that frustration freeness implies that the evolution should be localized “around the boundary”, and 
that for short times nothing at all would happen in the “bulk” (where the state is left invariant by the local 
interaction terms, because of frustration freeness). This intuition is formalized in the following lemma 

Lemma 12. Let C = (A4, B) an uniform family of Lindbladians, satisfying frustration freeness. Let A c T be 
a finite region, and fix a positive natural number m. Let B = A(m + 1), R = A(m + 1) \ A(m) and p ™ a fixed 
point ofT t A(m) and r an arbitrary state on R (see figure 1). 

II(if - T^ A ) (pZ ® < poly(m)W 1 (m) [e* - 1 + t]; (13) 

where T t BNj4 denotes the evolution generated by 

Z M Z+ Z B d d B . 

Z<zB\A d<m +1 


Figure 1: The construction of the sets B and R. 
B 



In order to prove such result, we will first prove a Lieb-Robinson-type of lemma. Denote by pit) = 
<g> r). For each X c B, denote L x the algebra generated by {Mz\Z c X}, which is the set of 
interactions terms of C B whose support is contained in X. 

Lemma 13. Under the same assumptions of Lemma 12, for each K e !Lx> the following “Lieb-Robinson-like” 
bound holds for some positive v: 

\\Kipitm, < poly(m)|X| \\K\l (e vt +1 - l)^ _1 (dist(X, R)). (14) 

Proof. Denote C(Z , t) - sup Te]Lz . Frustration-freeness implies that C(Z, 0) is 0 if Z n R = 0 (since 

Tr#p(0) = p ™), while is bounded by 1 otherwise. Moreover, let A(r) = ffd>r ||6^ s || o and for each Z c B, let 
6(Z) = A(dist(Z, R)). Assumption (A-2) implies that sup r v{r) A(r) < poly(?n). 

We are now going to replicate the proof technique of Lieb-Robinson bounds: denote C x = C B - C B '' X , 
and notice that, since they have disjoint support, [K, £ BnA ] = 0. Then 

= KC S p{t) = C B '' X K{p{t)) + KC x {p{t)), 

and consequently, 

t 

K(p(t))=e tcB ' X K(p( 0)) + J e {t - s)cB " X I\C X p{s) ds. 

0 

By taking norms 

t 

WXipim, < ||^(p(0))|| 1 + \Kl f ||£ x P ( S )| 1 d S ; 

0 

and thus 

t t 

CiX,t)<CiX, 0)+ y f \\MzlCiZ,s)ds + Y, [ K B || 0 ds 

ZrX±0 q d>m +1 q 

t 

<CiX,0)+5iX)t+ y \\Mzl f CiZ,s)ds. (15) 

ZnX±0 q 

By recursively applying equation (15), we obtain that 

oo j.n 

CiX,t)<Y J [a n + b n ]- 

n=n nl 
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where a 0 = C(X , 0) 


On= E ••• E \M Zl \\..4MzJC(Z n ,0), 

Z^X^O Z„nZ„_i*0 

6o = 0, b\ = 5(A) and 

b n+ 1 = E ••• E \\M Zl \\...\\MzAKZ n ). 

ZinXtO Z n nZ n -i±0 

Let us bound the two coefficients independently. The coefficients a n are threated in the same way as is done 
in the standard proof of Lieb-Robinson bounds [18]: recalling that C(Z , 0) is zero unless Z n R t 0, we have 
that 

a « = E E \M Zl \...\Mz n \ 

Z-^nX^O Z n rZ n - 1*0 
Z n rR±0 

We have that a\ is bounded by using assumption (A-l): 

°i - E E \\M Zl \\< 

ieX Z±3i 
Z\C\R±Q 

E W 1 (dist(i,R)) E lA^Zil J / (diam(Zi)) 

ieX Z\3i 

ZinRtO 

< v E W x (dist (i,R)). 

itX 


Similarly, we bound 02 as follows: 

®2 ^ E E E E ll-^Zi II I|A1 z 2 II • 

ieXZi3ijeZ2 Z^j 
ZrYtO 

We now use the fact that ^(dist(«, j))^(dist(j, R)) > jz(distand thus 

a 2 < E ^ _1 ( d ist(i,i?)) E \\MzA E Kdist(i,j)) E ^( dist (j,^)) \M Z *\ 

ieX Z\3i jzZ 2 Z 2 ^j 

ZnYt 0 

< E ^ _1 (dist(i,ii)) Ei Il-A^zil rz(dia m (Zi)) E E v(dieim(Z 2 ))\Mz 2 \ 

ieX Z\3i jzZi Z 2 ^j 

ZnYt 0 

< v E W^dist^,#)) E ||A^Zi|| rz(diam(Zi))|Zi| 

ieX Z\3i 

< v 2 E i /_1 (dist(i,i?.)). 

itX 

Proceeding in a similar way, we can bound a n by v n Y,izx v _1 (dist(i, R)). Let us now turn our attention 
to b n . Let Z\ n X + 0. Then for all u in Z\ it holds that 

dist(A, R) < dist(u, X) +dist(w,i?) < diamZi +dist(w,l?). 

In particular, this holds for yz 1 e Z\ such that dist(yz 1 ,i?) = dist(Ai, R). Therefore we have that 

1 < ^(diamZi)^(dist(Zi,i?))i/ _1 (dist(A',i?)). 

We can use the previous inequality to bound 6 2 as follow: 


6 2 = E WzAKZi) 

Zin.Y^0 

< iz _1 (dist(A, R)) E E P^Zil i / (diamZi)A(dist(Zi,i?))zz(dist(Zi,i?)) 

xeX Zi 3 x 

< upoly(m)|A'|r/^ 1 (dist(A, R)). 
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For 63 , we reason similarly as follows: 

&3= E I \\M Zl \\\\M Z2 \\5(Z 2 ) 

Z in.Y*0 Z 2 nZi*0 

IlMzJI 


I 


E E ||Mz 2 || i/(diamZ 2 )A(dist(Z2,i?))^(dist(Z2,i?)) 


0inJf*0 ^(dist(Zi,i?)) z z Zl Z 2 3 

<vpoly( to) E ||MzJ W 1 (dist(A L ,i?))|X L | 

ZinX *0 

< vpoly(m)^ _ 1 (dist(X,i?)) E E ||M Zl | |Zi|iy(diamZ 1 ) 

rreX Zi9a; 

< u 2 poly( 77 i)i/ _ 1 (dist(A, i?))|X|. 

Following the same argument, we can thus bound the general term 6 n+1 by 

v" poly(m)|X|j/ _ 1 (dist(X, R )). 

We can then bound 

t n 

Y a n~, < ( e vt - 1) E ^ 1 (diBt(*,iZ)) < (e vt - l)|AT|^- 1 (dist(AT, J2)) 


n\ 


ieX 


and 


t 


E^n—r^poly (m)v 1 (e vt -l)\X\u 1 (dist(X,i?)) + f[<5(A') - poly(m)|X|i/ 1 (dist(X, i?))]. 


Note that, because of Assumption (A-2), the last term in the r.h.s. can be bounded as 
(5(A) - poly(TO)|X|i/ - 1 (dist(A', R)) < poly(m)|X|i/ _ 1 (dist(X, R)) 
This concludes the proof. 

We can now prove equation (13). 

Proof. Applying Duhamel’s formula [40] we have that 

t 

( T B _ T B.A} {fC 9 t) = J T B^A£A p{s) ds . 

0 

t 

\\(T t Z-T t ^ A )(pZ®T)\\ < E f \\M z p(s)\\ 1 ds+ E 

1 ZnA*0f 


D 


and therefore: 


II B 9 d B \\t. 


d>m +1 


The second term on the r.h.s. is bounded by tpoly(m)^ ^m). Let us focus on the first term on the r.h.s. If 
Z c A(m), we can bound the r.h.s with equation (14). In particular, we have the following 

E ||AIzp(s )|| 1 <poly(ro)(e ws +t —1) E ll-^IL ]Z|i/ _ 1 (dist(Z, R)). 

ZnA*0 0nA*0 

ZcA(m) ZcA(m) 

Observe that, since Z n A + 0 , it holds that dist(Z, R) + diam Z > m, and therefore 

E ||A4z || 0 |Z|j/ _ 1 (dist(Z,i?)) < z/ - 1 (m) E ||Alz || 0 \Z\v{dvamZ) < , 

ZnA *0 0nA*0 

ZcA(m) ZcA(m) 

where we have used the Lieb-Robinson assumption. 
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If Z A(m), then it must hold that Z n R + 0. Then we showed in the previous lemma that: 


X ^ v\A\v 1 (dist(^4, R)) = v\A\v 1 (m). 

ZnAtfd 

ZnR±e> 


Putting it all together, we have that 

| (If - If’ ' A ) (pZ ® r)\l < poly(m)z,- 1 (m) [e* - 1 + t]. 


□ 


4 Decay of correlations 

In this section we show that, as a straightforward consequence of Lemma 11, the hypotheses on C imply that 
its fixed points have a particular character: they have fast decay of correlations, meaning that the correlations 
between two spatially separated regions is fast-decaying in distance. How fast this decay is is given by the 
decaying function A 0 defined in Lemma 11. 

Theorem 14. Under the same assumption as in Lemma 11, fix two regions A and B c A, let d A B > 0 be the 
distance between them. Then we have that 

T(A:B)<3(\A\ + \B\) s A 0 ^y (16) 

where the correlations are calculated with respect to p 0ol and 5 and Ao are defined in Lemma 11. 

Proof. Let C = Au B, and denote by pab the reduced density matrix of p a0 over C, and by pa and ps the 
reduced state on A and B , respectively. 

Consider p^ the unique fix point of T^ s \ and denote by p s A and p s B its reduced density matrices over A 
and B , respectively. If s < then C(s) has two disjoint components corresponding to >l(s) and B(s), 
and thus p decomposes as a tensor product over such bipartition, and its reduced density matrix over C is 
given by p s A <S> p s B . 

By Lemma 11, we have that, for any observable Oc with operator norm equal to 1 and supported on C, 

|Tr O c (p°o - pDI ^ |C| d A 0 (s). 


This implies that 

Wpab ~ Pa ® Pb\i - |CfAo(s). 

Since the trace norm does not increases under the partial trace, then 

lpA-p a Ah<\C\ S Ao(8), 

and the same holds for B. This in turn implies that 

II Pa ® Pb ~ Pa® PbWi - 2|C , |' 5 A 0 (s), 

and by applying the triangle inequality we obtain the desired result. □ 

Observation 15. By equation (4), we have that the mutual information I(A : B ) decays with dist(H,H) at 
essentially the same rate as T(A : B). 
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Note that the dependence on ]A| and \B\ of the bound in equation (16), which is harmless when A and B 
are finite regions - as in the case of two-point correlation functions - becomes significant when one of the two 
regions is proportional to the system size. In [44] the authors defined a correlation length £ for the mutual 
information, as the minimal length such that, for all L > £, it holds that 

I(b x (R - L) : b x (R) c ) < l -I{b x (R) : b x (R) c ), V®, VI? > 0. 

From this property they are able to derive an area law for the mutual information of the type: 

I (A : A c ) < 4|<9A|£. 

If we were to use (16) to determine a correlation length £, we would obtain a value for £ which depends on 
the system size, and thus obtain a bound comparable to equation (1) obtained in [20]. In the next section we 
want to show that it is possible to greatly improve this bound, at the cost of adding some extra hypotheses 
on the evolution. 


5 Area law for mutual information 

Until now, we have avoided as much as possible to put restrictions on the function v(r) appearing in 
Assumptions (A-l) and (A-2), and we have simply required it to be fast enough for Lieb-Robinson bounds 
to hold. Indeed, in most of the bounds obtained in the previous results it appears v(r) (or Ao(r), which 
depends on j/(r)), so that weaker assumptions will simply lead to weaker bounds. 

This gets more complicated when it gets to prove the area law bounds, since we are actually interested in 
pinning down the case in which the bound on mutual information takes the form of \dA\ log] A| - a area law 
with a logarithmic correction. As will be clear in the proof, the logarithmic correction depends on v(r ) (and 
consequently Ao(r)) to be exponential. Slower rates will still lead to a bound on the mutual information, 
but where the logarithmic correction is replaced by a super-logarithmic correction. The resulting bound can 
hardly be called an area law - a sub-volume law would be more correct. We will avoid such generalization, 
as they make the proof unnecessarily complicated, and focus on the more interesting case of exponentially 
decaying interactions, for which we can prove a proper area law with logarithmic correction. 


5.1 The pure fixed point case 

In this section, let us suppose that each T A has a unique fixed point, and moreover that this fixed point is 
pure. This setting is of particular interest in view of dissipative state engineering [23, 41], since ideally one 
would like to be able to create pure states (at least before noise and errors are taken into account). 

Let us denote by \(/) A ) the pure fixed point of T t A . 


Proposition 16. Let C = {A i,B} be a uniform family of dissipative evolutions, satisfying rapid mixing and 
having a unique pure fixed point \<j)\) for each A. Fix Ac A, and let pa denote the reduced density matrix of 
\cj) A ) on A. Then it holds that 

S(pa) < clog|A| • \dA\ 


for some constant c > 0. 


Proof. Let i > 0, to be determined later, and denote by a the reduced density matrix of \<t>A(e)) on A. Then 
we have trivially that 

S(a) <log 2 dim H A (e)\A ^ c 0 e\dA\, 

for some positive constant cq. On the other hand, by lemma 11, we have that d := \\p A - ct]^ < |A| < 5 Ao(I), and 
thus by equation ( 2 ) 


13 


S(pa) < S(cr) + iS'(cr) - S , (p j4 )| < 

Co t\dA\ + 2<i|^4| + 2h b (d) 

< c 0 £\dA\ + 2\A\ S+1 A 0 (£) + 2h b (d). 

Fix a e such that 0 < e < 1/2 and let us choose £ such that 2\A\ S+1 Aq(£) < e. This implies that £ scales as 
log|A|, and thus S{pa ) < C\ log|A| • \dA\ + e + 2hb(e/2\A\), for some positive Ci. By taking c > c\ we can absorb 
the terms depending on e in the other one, and obtain the claimed estimate S(pa ) < clog|A| • \dA\. □ 

While of interest, the case of a pure fixed point is a very specific one. Therefore, we want to give results 
applicable in the generic case of a mixed fixed point. To obtain such results we will need to make an additional 
assumption, namely that the system is frustration free. 

5.2 The frustration-free case 

Theorem 17 (Area law for mutual information). Let C - {At, B} be a uniform family of dissipative evolutions, 
satisfying rapid mixing, frustration freeness and having a unique pure fixed point. Let poa be the fixed point of 
C A for some A = b u (L) c T. Then we have that 

I(A:A c ) p „<c\dA\log\A\ (17) 

for some positive c independent of the system-size. 

Proof. For each n > 0, let p%^ be the fixed point of T A( ' n ' > . Fix a positive no to be determined later. Then it 
holds that 

I(A:A c ) p „ = 

I (A : A c ) p ~0 + X [I {A : A C ) PZ+ 1 - I (A : A c ) p „ ] . (18) 

n=no 

We want to show that it is possible to choose no in such a way that I (A : A c ) p ™ o < c|ch4|log|A| and the sum 
in the r.h.s. is arbitrarily small. 

For each n > 0, we have that, by rapid mixing (11) 

|| p n + i_ T A(„ + i) (p „ 0t) ||^ < 

where 4>\{n ) is a polynomial in n. On the other hand, equation (13) implies that 

\\T t Mn+1 ) (pZ ® r) - T t A ^ A (pl ® t)1 < |A|0 2 (n)W 1 (n)[e^ +1 - 1 ], 
and <j> 2 (n) is polynomial in n. Let us choose choose t n such that 

e n := </>i(n)e _7t " + ^(n)^" 1 (n)[e vtn + t n - 1] 

is exponentially decaying in n. This can be done by taking t n which scales proportionally to log ^n(n) 
which is essentially linear if is(n) grows exponentially. We can put the two bounds together using the triangle 
inequality, in such a way that 

|pr 1 -< (ri+ 1 )vA (^®T)| i <|A| e „. 

Observe that, since jA( n+ 1 A - 4 d oes no t ac t on A, /(A : d c ) r i(„+i)vA^„ ^ < I(A : A c ) p 
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Let us assume that no is big enough so that e n \A\ 2 <2 n for all n > no, we can then apply inequality (4) 
and obtain 


I(A:A c ) p ^-I(A:A c ) pZ < 


< I (A ■■ A C ) pS +l - I(A : A C ) T A(n + l)sA^ z ^ < 


6£n|^4|“ + Ah b (e n \A\). 


Then 

I {A : A c ) p „ < I(A : A c ) pZ0 + 6 Z 2- n + 4 Z h b (e n \A\). 

n=no n=n o 

Observe that, if 0 < x < 1/e, then ( x - 1) log 2 (l - x) < -a: log 2 x, and the latter is a increasing function in 
that interval. Therefore: 


h b (e n \A\) < -2£ n \A\\og 2 (e n \A\) < 


2~n+l 

-j^j-- + log 2 |^l|) . 


Therefore, Y,n=n 0 h b (e n \A\) is the tail of a series which is converging geometrically, and therefore is exponentially 
decaying as n 0 increases. The same is true for Y,n=n ^ n > so that both of them can be made smaller than 

i(a-.ac) pZ o. 

By taking no proportional to log|A|, we can bound I (A : A c ) p ™o by the logarithm of dim which is 
proportional to \A{uq) \ A\, and therefore: 


I(A:A c ) p2 o<c\dA\\og\A\. 


In conclusion, we have bounded the r.h.s. of (18) by c|ck4] log|A|, and this concludes the proof. 


□ 
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